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Abstract 

We study the realization of cosmic inflation in bigravity theories. By analyzing the 
evolution of scalar, vector and tensor perturbations in de Sitter-like spacetimes, we find 
strong stability constraints on the class of viable vacua offered by these theories. More 
specifically, the only stable de Sitter vacua contain two nondecoupled gravitons (one of 
which is massive) with different maximal propagation speeds. We derive an effective theory 
for the massless graviton, which is found to propagate at an intermediate speed, limited by 
the two maximal values. For inflation, while the spectrum of density perturbations remains 
nearly scale invariant, the power spectrum of tensor modes is found to depart from the usual 
prediction found in standard slow-roll inflation. In particular, both the tensor to scalar ratio 
r and the spectral index of tensor modes nx receive sizable contributions from the couplings 
of the theory, leading to specific signals that may be tested in future cosmological probes of 
CMB-polarization. 



1 Introduction 



Our present view of the universe relies entirely on the validity of Einstein's General Relativity 
(GR), which so far continues to be our best explanation to all known gravitational phenomena. 
Our confidence in GR is such, that from it we infer the existence of dark matter and dark en- 
ergy, required to fit various astronomical and cosmological observations, such as galaxy rotation 
curves, supernovae redshift-distance relation, CMB and large scale structure, just to mention a 
few. Despite of its undisputed success, a lack of a deeper insight on the nature of these dark 
substances still rises the question as to whether GR constitutes the correct theoretical framework 
by which the gravitational interaction should be addressed. 

Instead, GR might turn out to be an effective description of gravity valid at intermediate scales 
that needs to be completed at both UV and IR scales, leading to a more fundamental theory pQ. 
One concrete possibility is put forward by bigravity [2], which asserts the existence of a second 
spin-2 particle -in addition to the usual graviton- nontrivially modifying the long-range action 
of gravity. Recent developments have made clear that bigravity models imply significant but 
consistent departures from GR at long and short wavelengths [SHU [5]. As a consequence, these 
models offer an alternative view of phenomena such as dark matter and dark energy, leading to 
interesting prospects for future tests on gravity 0[7]- 

A key test on bigravity is whether it is able to incorporate an explanation to the origin of 
primordial density fluctuations [8] as observed in CMB [9\ and large scale structure measure- 
ments [10]. The purpose of this note is to address this question by analyzing the realization of 
cosmic inflation in bigravity models. We do this by first studying the evolution of perturbations 
on the class of de Sitter backgrounds offered by these theories. We show that only a restricted 
family of de Sitter vacua is stable under perturbations, characterized by the fact that its two 
gravitons remain coupled at wavelengths comparable to the horizon. By assuming that the de 
Sitter geometry evolves adiabatically towards a Minkowski vacuum (a quasi-de Sitter state) , we 
find that the power spectrum of tensor modes receives contributions that makes it differ from 
the usual prediction encountered in slow-roll models of inflation, offering a unique opportunity 
to test bigravity, involving measurements of CMB-polarization. 



We begin our discussion by introducing the basic setup to be studied, namely a system consisting 
of two metric fields g^ v and q^ v with inverse fields g^ u and q 1 *" respectively. The action describing 
this system is given by 
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where R g and R q are the Ricci scalars constructed from g^ u and q^ v , and X g and X q play the role of 
cosmological constants for each sector. The lengthscale L is introduced to make all parameters 
dimensionless, and may be taken as the fundamental lengthscale of gravity. The constant /3 
couples g and q with the help of a mixed volume element d 4 x(—q) u (—g) v with it + v = 1/2. The 
presence of /3 breaks the diff 2 -symmetry of the noninteracting theory down to diff [TT], making 
dU) to be only invariant under simultaneous gauge transformations of both metrics. We add to 
action (pQ) an interaction term first proposed in [2] 
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where k is the interaction strength. To keep our discussion simple, we consider the case u = 1/2, 
but point out that our results remain unchanged for the choice u = 0. Then, the Einstein's 
equations derived by varying flTJ) and ((2J) with respect to g pu and q pv are respectively given by 
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where G^ u (g) and G fMU (q) denote Einstein's tensors for each sector. Additionally, Tj{ v and T^ u 
are to be understood as the stress energy tensors sourced by g and q respectively when both (5 
and k are nonvanishing. These are given by 
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Finally, T^ 9 ^ and T^' represent stress energy tensors from matter fields coupled to each 
sector. 

We are interested in studying homogeneous and isotropic backgrounds. Disregarding intrinsic 
curvature effects, we may choose our two metrics to satisfy the following ansatze consistent with 
these requirements 
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a 2 (-dT 2 + dx 2 ), 
-X 2 dr 2 + y 2 dx 2 , 
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where a = a(r), Y = Y(t) and X = X(r) are scale factors that depend only on r, the conformal 
time with respect to g^. Notice that by a suitable change of coordinates, we can always find 
a frame where q pv is conformally flat instead of g^. For most of this discussion we focus our 
attention on vacuum solutions whereby T, 
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n{mq) 



0, and comment on the inclusion of 
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matter fields later on. We can already learn much from (|7|) and ([8]) by independently combining 
the 00 and 11 components of ([3]) and @ respectively, leading to 

a,; n,2 (X 2 -Y 2 )(/?y 2 -W) 

H' Y - HyHx = ~a 2{ - (10) 

where % a = a' /a, H x = X'/X, and %y = Y'/Y (here primes ' denote derivatives with respect 
to r). Now, the only way of achieving vacua characterized by H' a — % 2 = H Y — HyUx = 
is either by having X = Y or Y = 2\J n//3a. Examples of such vacua are precisely Minkowski 
and de Sitter vacua. The first branch X = Y corresponds to a case where both metrics are 
conformal to each other {q^ p oc gnu), whereas the second branch exists only if both k and f3 
are non- vanishing. To further understand these two branches we try the following scale factors 
representing two copies of de Sitter spacetimes 

a = — X = C x cl(t), Y = C y cl(t), (11) 

(with r < 0) where H$ is a positive dimensionless constant determining the expansion rate 
H = Hq/L of the homogeneous spacetime. Here, Cx and Cy are positive constants determining 
the second metric in terms of a{r). Given (fTT|) . eqs. ([9]) and (flO|) are simultaneously reduced 
to (C\ — C Y )(C Y (3 — 4k) = 0. There are only two additional independent equations, given by: 

3C X H$ + 6C y k - C x X g - C Y f3 = 0, (12) 
(3C| + C Y )P - 4C y (3tf 2 - \ q C 2 x ) = 0. (13) 

Let us start by analyzing the first branch Cx = Cy, referred to as the proportional vacuum 
(PV). In this case, it is straightforward to solve (|12p and ()13|) to obtain 



where we have defined X g = X g — 6k. Notice that Minkowski spacetimes may be obtained 
from this solution by tuning the parameters to satisfy X q X g = /3 2 . Additionally, in the limit 
K, /3 — > both metrics decouple and the solution reduces to standard de Sitter spacetimes for 
each gravitational sector. The second branch is called the non-proportional vacuum (NPV) and 
satisfies Cy = 2a/k//3. It exists only if k and j3 are non- vanishing, which implies that this 
solution cannot be reduced to a conventional de Sitter background by continuous deformations 
of the couplings. We find it convenient to express Cx and Hq in terms of a single parameter 6 
as: 

Cx = 9Cy, Cy = 2^, 2>Hl = X g -2^, (15) 
where 6 may be obtained by solving the following cubic polynomial equation: 

(3/3 + 16A g K//3) 6 3 + (/3 - 4X g )8 + 8k = 0. (16) 
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Observe that in order to have a Minkowski vacuum Hq = 0, one requires 9 = 2k/X 9 , and the 
previous equation implies that the parameters of the theory must satisfy (12k 2 +A 2 )/3 2 +64k 3 A, 2 = 
0. A remarkable property of the NPV is that the light cones of both metrics do not necessarily 
coincide. This means that each metric implies different maximal speeds to signals following 
their geodesies. While the maximal speed for particles following (/-geodesies is normalized to be 
c g = 1, the maximal speed for particles following g-geodesics is c q = 6. (It was shown in [I] that 
causality is preserved under these circumstances). Notice, however, that we could have worked 
in a frame where c q = 1 and c g = 1/0. 



3 Scalar and vector perturbations 

We now study the evolution of perturbations on these backgrounds. We proceed by expanding 
the two metrics as g^ = g^ v + h^ u + + h*, and g M „ = + r J, + r v v + r£, where g°„ 
and q® u represent the background solutions of eq. (fTTI) . and S, V and T denote scalar, vector 
and tensor modes respectively. Let us start by analyzing scalar perturbations. Since the theory 
is invariant under simultaneous gauge transformations on both sectors, we are then allowed to 
choose Newton's gauge to simplify one of the metrics, say g^ u , in the following way 

ftgj = -2a 2 ^, hfj = -2a 2 ^5ij. (17) 

Since this gauge leaves no residual symmetry, r^L, cannot be reduced in any similar way, and 
must be treated in its most general form: Tq = — 2X 2 A, rj| = rf = XYdiF, and rf- = 
Y 2 (-2BS ij +d i d j E), where A, B, F and E are the four scalar modes of the g-metric. By inserting 
these perturbations back into (j3J) and (j4|) we deduce the linear equations for the evolution of all 
scalar modes. In the particular case of the PV solution Cx = Cy, the combination £ = + ^ 
decouples from the rest of the modes and satisfies the following equation of motion 

£" + -S' + A: 2 S + M 2 (r)S = 0, (18) 
r 

where k labels the modes wavenumber in Fourier space. The mass M 2 (t) in eq. (|18p is found 
to be: 

M 2 _ /3 + X 2 2((3-X 9 + X q ) 2((3-X g ) 

2a? 3L 2 x ^ 3L 2 ^ 3L 2 C 2 X ' [ 1 

Because of the sign of the second term in (|18p . the combination $ + grows without control 
regardless of the values of /3, k, X g and X q , rendering PV solutions unviable to accommodate 
inflation. This result was already found in [7] for the particular case k = 0. There, it was 
speculated that with the inclusion of Sint it would be possible to cure this instability. As we 
have seen this instability persists, however, the existence of NPV solutions for k / opens up 
the possibility of having new stable backgrounds for the propagation of perturbations. This is 
indeed what we find: in the case of the NPV solution (|15p all scalar perturbations are constrained 
to vanish (satisfying equations of the form fc 2< 3> = 0) meaning that this background is purely 
geometric. This is corroborated by the fact that vector perturbations and rY v also vanish 
in this vacuum. 
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4 Tensor perturbations 



The previous discussion implies that the only viable vacuum to realize inflation is the NPV 
solution, where the only dynamical degrees of freedom are the two tensor modes hj- and rfj, 
which we now study. To proceed we write hjj = a 2 ^ and rf- = a 2 V9CyXijj with 7^ and 
Xij traceless and restricted to satisfy <9 l 7ij = d l Xij = 0. The factor \fdCy in front of Xij has 
been introduced to ensure both fields 7^ and Xij to have the same kinetic energy normalization. 
Disregarding indices, the equations of motion for 7^ and Xij are found to be: 



„ 2 2 C 2 a 2 f3(9 2 - 
7 '-1+^- — jt 

n 2 , o2l2 a 2 f3(6 2 - 
X - -X + G 2 k 2 X + H \ 2 
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(20) 
(21) 



This set of equations describes a coupled system of gravitons, of which, only one is massive. The 
nonzero eigenvalue of the mass matrix is given by 



a 2 p(e 2 -1)(0 + C Y )/9L 2 , 



(22) 



from where we obtain the restriction 9 > 1 if (3 > and 9 < 1 otherwise. Observe that in the 
short wavelength limit rk 3> 1, g-gravitons propagate with speed c g = 1, whereas g-gravitons 
propagate at c„ = 9. However, for k < m a both modes remain mixed, and we have to proceed 
carefully. The massless and massive modes, hereby denoted £ and <r, may be obtained from 7 
and x through a rotation as: 



v^7 + Cyx 
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+ C 2 
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The equations of motion for £ and a are then given by: 



f " - + A; 2 %£ + fc 2 ^ = 0, 

2 

<x" - -o-' + k 2 A aa a + A; 2 ^£ + m 2 a a = 0, 



(24) 
(25) 



where % = 0(1+C 2 0)/(fl+C 2 ), ^ = (03 +C 2 )/(0+C 2 ) and = = V9{9 2 -1)C Y /(9+ 
Cy). The two fields continue to be coupled through their kinetic terms. However, if m 2 (r) S> 
2/r 2 , the massive graviton decays quickly before horizon exit, and the only relevant degree of 
freedom becomes the massless mode. One way of addressing this situation is by deducing an 



effective theory for the massless mode valid for the regime k <C mi.. Following |12j . we find 



that the massive field may be expressed in terms of £ as a : 
and the dynamics of the massless mode is well described by 



-k 2 A a £/(m 2 a - 2/r 2 + k 2 A aa ), 



-t' + k 2 
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2/r 2 + k 2 A 



e = o. 
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This corresponds to a massless graviton with a modified dispersion relation. Notice that if the 
parameters of the theory X g , X q , /3 and k are all of order unity, then the condition m 2 3> 2/r 2 
is equivalent to Hq <C 1 (or H <C L^ 1 ), which is necessary in order to trust the present field 
theoretical description of our system. Since we are interested in phenomena for which k 2 <C m 2 
the contribution coming from the second term inside the bracket in (|26h may be neglected, and 
the speed of propagation Ch for this mode becomes 

cl = e(i + c 2 Y 0)/(e + c 2 Y ). (27) 

Notice that Ch G [1,0], depending on the value of C Y = 4k/ (3. This is because the massless 
mode appears from the mixing between both metrics, and therefore its propagation is affected 
by both backgrounds simultaneously. 



5 Discussion 

The previous results have some interesting and nontrivial consequences that we now discuss. To 
achieve realistic models of inflation we need to move our analysis from de Sitter to quasi-de Sitter 
spacetimes, in such a way that the background quantities Hq, Cx and Cy evolve adiabatically 
towards Minkowski. This may be achieved, for instance, by introducing a scalar field <j) into the 
theory, and letting A 9 , X q , (3 and k to be functions of it. The dependence of these parameters 
on (j> will be restricted to satisfy certain slow-roll conditions, just as in the case of conventional 
slow-roll inflation. Disregarding these model dependent aspects, since in pure de Sitter there 
are no scalar degrees of freedom, in this new quasi-de Sitter phase there will necessarily exists a 
curvature mode £ reflecting the time translation symmetry breaking, which under very general 
assumptions [13] will be described by the following effective quadratic action: 

S( = ~~hj d 3 xdT^e[((T-c 2 s (d() 2 ], (28) 

where e = —H/H 2 and c s is the speed of sound for density perturbations, which is expected 
to depend on Cx and Cy (here ' denotes a derivative with respect to cosmic time dt = adr). 
This situation is rather generic in the context of inflationary physics, and leads to the well 
known power spectrum for scalar densities: Vs{k) = L 2 H 2 / (8ir 2 ec s )\k Ca =aH , with spectral index 
n s = 1 + r) — 2e — e s where r\ = —e/He and e s = c s /Hc s . 

The novel aspect comes with tensor modes: since they propagate with a variable speed Ch 
in the regime c^k 2 <C m 2 , the power spectrum of tensor modes departs from the conventional 
prediction encountered in GR. To arrive to a concrete prediction, we assume that matter fields 
are coupled to and in such a way that it only couples to the massless mode £ [5]. 
Other alternatives may be considered but would lead to similar conclusions. Then, by imposing 
Bunch-Davis vacua for subhorizon modes (in the regime a 2 H 2 <C k 2 -C m 2 ) we find: 



V T(k) = -jrg- 



(29) 

kci,=aH 
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This result modifies the expression for the tensor to scalar ratio to be r = 16 e c s /c? . Additionally, 
the spectral index tit of tensor perturbations is now found to be tit = — 2e — 3e/j, where = 
dh/Hch. Since is sensitive to the specific dependence of A 9 , X q , (3 and k on cj), the power 
spectrum Vt{^) m bigravity models is no longer restricted to be red tilted. Finally, the effects 
of the UV field a on 'Pr(fc) are of order H 2 /m 2 <C 1. 

Thus, we see that the realization of inflation severely constrains the parameter space of bigrav- 
ity models. Our main result is sumarized in eq. (j26|) which shows that tensor modes propagate 
with a modified dispersion relation. As a consequence, the predicted power spectrum for tensor 
modes generically differ from the conventional prediction offered by slow-roll inflation. In addi- 
tion, we expect other relevant departures from the conventional picture such as the enhancement 
of non-Gaussian distribution of tensor modes (see also |14j). The importance of these results 
is twofold. On the one hand, it forces us to widen our view of effective theories of inflation 
to include consequential modifications to the tensor sector of the theory, consistent with the 
symmetries of quasi-de Sitter spacetimes [13|. I15j . On the other, it re-emphasizes one more 
time the importance of improving the precision of CMB polarization measurements |16| to test 
physics of the very early universe. Since primordial B-mode signals are exclusively due to tensor 
perturbations, the next generation of CMB polarization probes might give a powerful insight 
into bigravity theories. Last but not least, notice that eq. (|26p is also valid for Minkowski space- 
times, and therefore it implies significant levels of departures from GR for gravitational waves 
phenomenology. For instance, c/j > 1 would involve a travel delay of gravitational waves from 
their sources when compared to light signals. Otherwise, Ch < 1, would produce a Cherenkov 
type of radiation emitted by particles exceeding c^. 
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